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Abstract We present a solution for the stationary state of an asymmetric ex- 
clusion model with sequential update and open boundary conditions. We solve the 
model exactly for random hopping in both directions by applying a matrix-product 
formalism which was recently used to solve the model with sublattice-parallel up- 
date 0. It is shown that the matrix-algebra describing the sequential update and 
sublattice-parallel update are identical and can be mapped onto the random sequen- 
tial case treated by Derrida et al@. 
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The one-dimensional asymmetric exclusion model (AEM) is one simple example of a 
reaction-diffusion model and has been used to describe various problems in different 
fields of interest, like the kinetics of bio polymerization Jl and traffic [|]]. Using 
recursion relations on the system size it was solved 1992 || for the case of random 
sequential update and open boundary conditions. Open boundaries here and in the 
following mean that particles are injected at one end of a chain of L sites with proba- 
bility a and removed at the other end with probability /?. This model was then solved 
again by Derrida et al 1993 |fj] using a matrix product Ansatz (MPA), inspired by the 
Matrix Product Groundstate for quantum spin chains with groundstate energy zero 
||, for the weights of the stationary configurations. This Ansatz was most elegantly 
used to obtain expressions for the density profile and higher correlations consisting 
only of products of two matrices (belonging to the two possible states of each site) 
and two vectors describing the influence of the boundaries. Since then the MPA was 
extended to find also the transient of the model J7| and to recover solutions of certain 
known integrable reaction-diffusion models ||. All these models work with random 
sequential update. 

Hinrichsen|| solved the AEM model with sublattice-parallel update and could con- 
firm earlier conjectures for the correlation functions ||. It should be made clear that 
this update is substantially different from the fully parallel update used for example 
for modelling traffic flow. Nevertheless, to our knowledge this is the first model with 
noncontinous time which has been solved using the MPA. It is an open question for 
which classes of models the MPA can be used successfully. 

In this letter we study the AEM with sequential update. The model then can be 
defined as follows: 

Consider sites located on a chain of length L. Each site i (0 < i < L) can be occupied 
by a particle (rj = 1) or may be empty (r,; = 0). We start the update at the right end 
of the chain and remove a particle at site i = L with probability (3. We then look at 
the pair of sites (i = L — 1, i = L). If we find a particle at site L — 1 and no particle 
(called a hole) at site L we move the particle one site to the right with probability 
p. In the opposite case the particle hops one site to the left with probability q. In 
the remaining two cases nothing happens. We continue the update with the pair 
(i = L — 2,i = L — 1) and so forth until we reach the left end of the chain. After the 
update of pair (i — l,i — 2) we inject a particle at site i = 1 with probability a if 
the site is empty. 

Let us remark at this point that this model put on a ring (no injection/removing 
of particles and periodic boundary conditions) has a trivial stationary state JTO, [TT| 
where correlations are absent. 

One could reverse the order of the sequential update and go from the left to the right 
through the chain (using the same rule as before for the pair-update). These two 
models are connected by a particle-hole symmetry: injecting particles can be seen as 
removing holes, and vice versa. Therefore it is sufficient to study just one model. The 
particle-hole symmetry gives a first hint at the phase diagram: one would naively 
expect the same phase diagram for both models. This implies a symmetry of the 
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phase diagram in a and (3. 

Following the MPA, we write the stationary probability distribution Po(tl, T2, Tl) 
as 



Po(n, t 2 , r L ) = Z L \W\ U(nD + (1 - r0^)| V) 



i=i 



or formally as 



|P ) = Z L -'(W\ (|j)° £ |V). 



(1) 



(2) 



The square matrices £ and £) can be infinite-dimensional. The vectors {W\ and |V) 
act in the same vectorspace as E and D and will, like E and D, depend on a and 
/3. The normalization constant is given by Zl = (W\(D + E) L \V). \Pq) represents 
the weights of all the configurations in the stationary state. This implies that \P ) is 
invariant under the action of the update-operator or transfer matrix T: 



TlPn 



(3) 



Let us now explicitly write down T. The boundary conditions may be represented 
by operators 1Z and C acting on site i — L and i — 1, respectively: 



K = 



1 p 
l-(3 



C = 



I- a 
a 1 



(4) 



The chosen basis for 1Z and C is (0, 1). The update-rule for any pair of sites (i, i + 
can be written as 



/ 1 
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The basis is (00, 01, 10, 11) and we have formally 

T — L - Ti • ... • • R 

with 



(5) 



L = C 
R = 1( 
Ti = lc 



> 1 <g> ... <g> 1 , 
... ® 1 ®K, 
1...®%® 1... 



1, 



(6) 



(7) 
(8) 
(9) 



where 1 denotes the identity matrix. 
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It is peculiar that one can use exactly the same mechanism which was used to deter- 
mine the stationary state of the model with sublattice-parallel update fll] in order to 
fulfil equation (3): 



r 



(EMI) =(D°(S) 



(10) 



(W\C 



\V) 



with some square matrices E, D. This means that a 'defect' is created in the be- 
ginning of an update at site i = L and then transported through the chain until it 
reaches the left end, where it disappears. 
Equation (10) leads to the following bulk algebra: 



[E, E] = [D, D] 
(1 - q)ED + pDE 
qED + (1 - p)DE 



and the boundary conditions 



(W\E(l-a) 
(W\(aE + D) 



(W\E 
{W\D 



, 
ED 
DE 



(1-{3)D\V) 
(E + (3D)\V) 



D\V) 
E\V) 



(12) 



For q = 0, p = 1 we recover the algebra which was solved in |l[] for the model with 
sublattice parallel-update, and we can adopt the two-dimensional representations of 
E, D, E, D, (W\, \ V) of the algebra (7), even though \P ) has a different structure. 
One can easily show that the density profile of the sequential update corresponds to 
the density of the even sites with sublattice-parallel update and one gets essentially 
the same phase diagram. 
One can check that for the case 



(13) 



one finds a one-dimensional solution of the complete algebra. This equation defines 
the lines in the phase diagram on which the mean field solution becomes exact. One 



can use these lines to calculate different currents in the phase diagram |10|, [TT| and 
we can exclude the case p = q = 1 in the following. 
To solve the general algebra, we first note that 



[E + D, E + D] = 
holds for all values of p, q. By demanding 



E 
D 



E + Xl. 
D — XI 



(14) 



(15) 
(16) 
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(with some real number A and the identity matrix 1) one can reduce the whole algebra 
of seven equations to just three equations: 



We define 



pDE - qED 

a(W\E 
0D\V) 



D 
E 

A 2 



A(l -q)E + A(l -p)D 

X(l-a)(W\, 

\\V). 



X{l-p)D 
X(l-q)E, 
1 



(17) 



18) 



(l-g)(l-p) 



and rewrite eq.(14) as 



pDE - qED 
a(l-p)(W\E 
P(l-q)D\V) 



E + D, 

(l-a){W\ 

\V). 



(19) 



This is the algebra for the AEM with random sequential update but the same local 
transfer matrix and the same boundary conditions as in our model. It was solved by 
Derrida et al with infinite-dimensional matrices. Note that one has to rescale the 
vectors (W\ and \V) of their solution with and yr^, respectively. 
For the case q = we write down a slightly different representation (A = 1): 



D 



1 

V 



( | ai 

110 

11 

1 



E 



1 

V 



( p(l-a) 
a 








1 — p 








p 1 — p 













p 1 — p 



\ 



(20) 



<W| = (1,0,0,0), \V) 



(I1CL2 



U {1 








a) (3 + a — p] . 



(21) 



(22) 



The MPA makes it now straightforward to calculate the density profile, higher cor- 
relations, the current and the phase diagram. We note at this point that for the 
case q = the randomness in p produces qualitatively the same phase diagram as 
in the case of random sequential dynamics with q — 0, p — 1. Further results, also 
concerning the comparison of the AEM with random sequential, sublattice-parallel, 
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sequential and fully parallel update will be published elsewhere |IOL ITT 



The present work shows that the sequential AEM with stochastic hopping in both 
directions can be solved by means of a mapping onto the known solution of the 
random sequential AEM. This implies also a generalization of the known solution 
of the AEM for sublattice-parallel update, which had been solved for deterministic 
hopping. It is remarkable that pDE — qED = E + D is the fundamental bulk al- 
gebra equation for all three types of update discussed so far. It is now possible to 
understand more about the nature and implications of different types of update and 
randomness [10 , |TT| . 



The solution of the sequential dynamics gives reason to hope that the MPA is capable 
of solving the fully parallel case, which can be written as a three-state model with 
the same sequential update as studied in this workJTJJ, [TT| . 
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